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(~| . Abstract. The Chandra LETG spectrum of SS Cyg in outburst shows broad 5 A) 

pi f spectral features that have been interpreted as a lar ^e number of a bsorption Unes on a 

' ' blackbody continuum with a temperature of 250 kK (Mauche 2004). It is most probable 

that this is the spectrum of the fast-rotating optically thick boundary layer on the white 
dwarf surface. Here we present the results of fitting this spectrum with high gravity 
hot stellar model atmospheres. An extended set of LTE model atmospheres with solar 
chemical composition was computed for this purpose. The best fit is obtained with the 
following parameters: Teff = 190kK,logg = 6.2,andA^H - 8-10''^ cm"^. The spectrum 
of this model describes the observed spectrum in the 60-125 A range reasonably well, 
but at shorter wavelengths the observed spectrum has much higher flux. The reasons 
for this are discussed. The derived low surface gravity supports the hypothesis of the 
fast rotating boundary layer. 



1. Introduction 

SS Cyg is one of the bes t-studied cataclysmic variable stars and is a prototype of dwarf 
nova stars ( Wameiiri995h . X-ray radiation of this close binary in quiescence is hard and 



can be described by an optically thin hot (IcT « 20 keV) plasma with an observed flux 
w 2- 10"'" erg s"' cm"^. In outburst, the hard X-ray flux decreases by a factor of ten, the 
plasma temperature is reduced to ~ 6-8 keV, and an additional soft component appears 
with a bla ckbodv temperature ~ 200-300 kK dCordova et all Il980l : iMcGowan et al. 
12004 : Jshida et alj|2009 ). 

A high-resolution spectrum of the soft c omponen t was o btained with the Chan- 
dra LETG and was carefully investigated by lMauch3 ( 20041) . Th e spectrum, which 



looks hke the pho tospheric spectra of super-soft X-ray sources ( Rauch et all l2010l : 



van RossumI |2012|) . is naturally associated with the radiati on of the boundary layer 



(BL) betw een the accretion disk and white dwarf (WD) dPringle & Savonije 1 19791 : 



Kleylll991l) . and can be phenomenologically described by a blackbody spectrum with 
T w 250 kK and numerous broad absorption features of ion s of O, N e, Mg, Si, S, and 
Fe; the BL luminosity and WD spin were also evaluated by lMauch3 (l2004h . Here we 
present our attempt to fit the Chandra LETG spectrum of SS Cyg using the spectra of 
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hot stellar model atmospheres that are close to the Eddington limit, and to make more 
accurate estimates of the BL parameters on this basis. 



2. Model atmospheres 

To model high temperature atmospheres th at are close to the Eddington limit, we used 
our version of the computer code ATLAS ('Kurucz"l970'), modified to deal with high 
temperatures (Ibragi mov et al. ,2003; S uleimanov &JWemer 2007 ) . We assumed local 
thermodynamic equilibrium (LTE) and ac counted for the pressure io nization effects 
using the occupation probability formalism ( Hummer & Mih alas ' 1 9 8 S") as described by 



iHubeny et al. (1994 ). We took into account coherent electron scattering together with 
the free-free a nd bound-free opacity of al l ions of the 15 most abundant elements using 
opacities from I Vemer & Yakovle 3 (11995). Line blanketing is taken into acc ount using 



25000 spectral hnes from the CHIANTI, Version 3.0, atomic database (iDere et al 



19971) . 

Using our code, we calculated 22 model atmospheres with solar chemical com- 
position. The effective temperatures of the models range between 150 kK and 250 
kK with a step of 10 kK. We used two values of the surface gravity for each effec- 
tive temperature: log^ = log^Edd + 0.2 and logg = loggedd + 0.4, where log^Edd = 
log(o"e cTsBr^jj/c) = 4.88 + 41og(reff/10^ ^) is the surface gravity that has an equal 
radiation pressure force for a given T^ff, and cr^ « 0.34 g cm"^ is the electron scattering 
opacity for the assumed solar chemical composition. The positions of the computed 
models on the Teff-logg plane are shown in Fig. [T] (left panel). The considered model 
atmospheres are very close to the Eddington limit and a radiation pressure force g^ad 
due to spectral lines becomes larger than the surface gravity at the upper atmospheric 
layers (see Fig.[T] right panel). We did not consider moving atmospheres and we took 
a gas pressure equal to 10% of the total pressure (Pgas - O.lPtot) if ^rad > ^ to enforce 
hydrostatic equilibrium. 

Examples of the computed emergent spectra together with the temperature struc- 
tures are shown in Fig.|2] The emergent spectra are dominated by absorption line forests 
and they become similar to the observed spectrum of SS Cyg only after accounting for 
the LETG spectral resolution (see Fig.[3ll. The differences between the spectra of mod- 
els with different surface gravities is obvious (see Fig.|2l right panel) and can be found 
from a comparison with the observed spectrum. 



3. Results 

We fit the observed soft X-ray spectrum of SS Cyg using the set of model atmosphere 
spectra described above, convolved with the Chandra LETG spectral resolution AA = 
0.05 A. The comparison of the best-fit spectrum with the observed spectrum is shown 
in Fig. 121 and the contours ofx^ on the Teff-log A^h parameter plane are shown in Fig.|4] 
The fitting procedure was performed in the 60-125 A wavelength range only, because 
at the shorter wavelengths our model spectra cannot describe the observed spectrum. 
The possible reasons for this are discussed in the next section. We found that the best- 
fit model parameters correspond to the models with the lower surface gravity log g = 
loggEdd + 0.2 with Teff = 190 kK, A^h - 8 • 10^^ cm^^, and the normalization K = 
fRlr^/d'^ ^ 7.82 • 10"^^, where d is the distance to SS Cyg and / is the WD fractional 
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Figure 1. Left: Computed model atmospheres on the Teff-log^ plane. The locus 
of points for which log g - log gEdd is shown by the dashed curve. Right: Ratio of 
the radiation force to the surface gravity vs. depth for various model atmospheres. 




Figure 2. Left: Emergent spectra and temperature structures of three model 
atmospheres with the same differences between \ogg and the Eddington limit: 
\ogg - log gEdd + 0.2, and various effective temperatures: 150 kK (solid curves), 
200 kK (dashed curves), and 250 kK (dotted curves). Right: Emergent spectra 
and temperature structures of two model atmospheres with the same effective tem- 
perature (200 kK) and different \ogg: \ogg = log gEdd + 0.2 (dashed curves) and 
log g = log gEdd + 0.4 (solid curves). 
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T(kK)=190, log(g) = 6.20, logNH=19.90, chisqr/dof= 3.92 

3.0 1 • • • \ if • • \ ^ • • \ • • • \ ' 




40 60 80 100 120 

wavelength (A) 



Figure 3. The Chandra LETG spectrum of SS Cyg in outburst (thick black curve) 
and the best-fit model atmosphere spectrum with Tefj =190 kK, logg = 6.2, and 
logA^H - 19.9 (thin red curve). The fitting was performed in the 60-125 A wave- 
length range. The model spectrum at the shorter wavelengths is shown by the dashed 
curve. 
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Figure 4. Position of the best-fit model on the Tetf-log A^h parameter plane and 
contours of^^ = [1.5, 3, 6] 
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T(kK)=190, log(g) = 6.20, logNH(cm"^)= 1 9.90, x/=7.65 
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Figure 5. The Chandra LETG spectrum of SS Cyg in outburst (thick black curve) 
and best-fit model atmosphere spectra with Tefr = 190 kK, logA^e = 19.9, and two 
surface gravities: logg - 6.2 (thin red curve, top panel) and \ogg - 6.4 (thin blue 
curve, bottom panel). The fitting was performed in the 45-125 A wavelength range. 



area occupied by the BL, which can be expressed as the relative BL extension along the 
WD surface / ~ (InRwD IHBiMi^nR^^) = Hbl/Rwd- 

The obtained fit is not completely satisfactory (reduced = 3.9), hence the formal 
errors are large and we have not attempted to determine them. However, the choice of 
the lower surface gravity is statistically justified because the reduced is significantly 
larger for fits obtained using the higher gravity models (Y^of ~ ^-^^ ^^-^dof ~ ^-^^ 
for fits performed in the 45-125 A wavelength range). This is also clear from Fig.[5j 
which shows that the high gravity spectrum cannot describe the local flux maxima at 
78 and 92 A. 

Using the obtained fit parameters, we evaluated some basic properties of the BL. 
For ease of co mparison, we adopted the same system and WD parameters used by 
Mauchd (|2004|) : Mwd = 1 Mq, /?wd - 5.5 • 10^ cm, d = 160 pc, and an accretion disk 



bolometric luminosity in outburst Loisk - 10^^ erg s ^ In this case, we can evaluate the 
fractional area of the BL f = 6.3 ■ 10"^ (5 ■ 10~^), the bolometric BL luminosity Lbl - 
1.8 • 10^"* (5 • 10^^) erg s-^ ar id the relative B L luminosity LeL/Z^Disk - 0.18(0.05), 
where the values obtained by iMauchd ( 20041) are shown in parenthes es. Using the 
well-know n relation between the BL and accretion disk luminosities (Kluz nial3ll987l : 
lKleylll99lh Z^BL/^oisk - [1 - ^^wd/^^k(^wd)]^» where Qk(^wd) is the Kepler angular 
velocity at the WD radius, we infer that the spin period of the WD in SS Cyg is 12 (9) s. 

The intrinsic surface gravity of the WD in SS Cyg (loggwD - 8.46) is more than 
two orders higher than the obtained BL effective surface gravity log geff - 6.2. The 
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surface gravity of the BL can be reduced by the fast rotation of the accreting matter, 
and we evaluated a BL angular velocity Obl * 0.98 Qk(^wd) using the simple relation 

geS = gWD - ^ujRwD - gWD (1 " [^Bl/^^k(^Wd)]^)- 



4. Conclusion and Discussion 

On the basis of the above analysis, we conclude that the BL in SS Cyg can be considered 
as a hot {» 190 kK), fast rotating [Qbl « 0.98 Qk(^wd)], narrow {Hbl * 0.063 /?wd) 
belt on the WD surface. 

This deduction is founded on the fit of the SS Cyg Chandra LETG spectrum 
with the model atmosphere spectra. The obtained fit is not statistically acceptable 
(jC'^^f = 3.9), especially at the shorter wavelengths (< 60 A) and in the 82-90 A wave- 
length region. These deficiencies can be connected with model shortcomings: e.g., the 
cher nical composition may differ from solar, non-LTE effects could be important (see, 
e.g.. iRauch et al.ll2010h . and the atomic data are almost certainly neither complete nor 



entirely accurate. The most important unmodeled effect is atmosphere expansion due 
to a spectral line driven stellar wind, which ca n be significant bec ause grad > ^ at the 
outer layers of our model atmospheres (see also Ivan Rossumll2012h . It is likely that the 



BL cannot be described by a simple one-zone model and that it has a more compli- 
cated structure, with a distribution of effective temperatures and surface gravities over 
its surface. All of these effects must be taken into account in further investigations. 
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